Anisotropic flow in Cu+Au collisions at y/s^N = 200 GeV 
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The anisotropic flow of charged hadrons in asymmetric Cu+Au collisions at the Relativistic Heavy 
Ion Collider is studied in a multi-phase transport model. Compared with previous results for sym- 
metric Au+Au collisions, charged hadrons produced around midrapidity in asymmetric collisions 
are found to have a stronger directed flow v\ and their elliptic flow V2 is also more sensitive to 
the parton scattering cross section. While higher-order flows V3 and V4 are small at all rapidities, 
both vi and vi in these collisions are appreciable and show an asymmetry in forward and backward 
rapidities. 

PACS numbers: 25.75.Ld, 24.10.Lx 



I. INTRODUCTION 

There have been extensive studies on the azimuthal 
anisotropy of hadron momentum distributions in the 
transverse plane perpendicular to the beam direction, 
particularly the elliptic flow vi, in heavy- ion collisions 
at various energies [if. The hadron transverse momen- 
tum anisotropy is generated by the pressure anisotropy 
in the initial compressed matter formed in non-central 
heavy ion collisions 0, Q and is sensitive to the proper- 
ties of produced matter in these collisions. For heavy- ion 
collisions at the Relativistic Heavy Ion Collider (RHIC), 
it has been shown that this sensitivity not only exists in 
the larger elliptic flow 0, 0, la m IE 111 lill liil but also in 
smaller higher-order anisotropic flows |l2tll3l Hil liH liH 
ITtI Il8l | . To investigate the influence due to initial colli- 
sion geometry on anisotropic flows in heavy ion collisions, 
one usually varies the impact parameter of a collision or 
the atomic number of colliding nuclei (l9|. Another pos- 
sibility is to study asymmetric heavy ion collisions with 
unequal mass nuclei as their initial collision geometry is 
different from that in collisions with equal mass nuclei. 

In the present work, we use a multi-phase trans- 
port (AMPT) model, that includes both initial partonic 
and final hadronic interactions |20l l2l| . to study the 
anisotropic flows v\, V2, V3, and U4 of charged hadrons 
in asymmetric Cu+Au collisions at y/s = 200 AGeV at 
RHIC. Use is made of both the default version and the 
version with string melting, i.e., allowing hadrons that 
are expected to be formed from initial strings to con- 
vert to their valence quarks and antiquarks 0> |2^, |23| . 
The latter was able to explain measured pt dependence 
of V2 and V4 of mid-rapidity charged hadrons with a 
parton scattering cross section of about 10 mb while 
the default version seemed to give a better description 
of the anisotropic flows at large rapidities where the 
string degree of freedom dominates |2J|. We find that 
charged hadrons produced around midrapidity in asym- 
metric Cu+Au collisions have a stronger directed flow v\ 
and their elliptic flow is also more sensitive to the par- 
ton scattering cross section than those found previously 



in symmetric Au+Au collisions. Furthermore, both v\ 
and i>2 of charged hadrons are asymmetric in the forward 
and backward rapidities. The higher-order flows V3 and 
Vi are generally small, and their sensitivity to the parton 
scattering cross section is thus less clear. 

This paper is organized as follows. In Sec. [HJ the 
AMPT model is briefly reviewed. Results on the pseudo- 
rapidity dependence of the anisotropic flows of charged 
hadrons are shown in Sec. IIIII The px dependence of 
the anisotropic flows of charged hadrons around midra- 
pidity is studied in Sec. IIVI while that at large rapidities 
is shown in Sec. Finally, a brief summary is given in 
Sec. EH 



II. THE AMPT MODEL 

The AMPT model [S3, E3, EE EE H3 is a hybrid model 
that uses minijet partons from hard processes and strings 
from soft processes in the Heavy Ion Jet Interaction Gen- 
erator (HIJING) model |2^ as the initial conditions for 
modeling heavy ion collisions at ultra-relativistic ener- 
gies. Time evolution of resulting minijet partons is de- 
scribed by Zhang's parton cascade (ZPC) [2<J model. At 
present, this model includes only parton-parton elastic 
scatterings with an in-medium cross section given by: 



dfj p 
~dt 
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where the strong coupling constant a s is taken to be 0.47, 
and s and t are usual Mandelstam variables. The effective 
screening mass [i depends on the temperature and den- 
sity of the partonic matter but is taken as a parameter in 
ZPC for fixing the magnitude and angular distribution of 
the parton scattering cross section. After minijet partons 
stop interacting, they are combined with their parent 
strings, as in the HIJING model with jet quenching, to 
fragment into hadrons using the Lund string fragmenta- 
tion model as implemented in the PYTHIA program [3(| • 
Final-state scatterings among hadrons are then modelled 
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by a relativistic transport (ART) model |3l| . The default 
AMPT model [2(| has been quite successful in describing 
measured rapidity distributions of charge particles, par- 
ticle to antiparticle ratios, and spectra of low transverse 
momentum pions and kaons | 2l| in heavy ion collisions at 
the Super Proton Synchrotron (SPS) and RHIC. It has 
also been useful in understanding the production of J / ip 
|25| and multistrange baryons |2q | in these collisions. 

Since the initial energy density in heavy ion collisions 
at RHIC is much larger than the critical energy den- 
sity at which the hadronic matter to quark-gluon plasma 
transition would occur [H IS El, the AMPT model 
has been extended to convert initial excited strings into 
partons |10| . In this string melting scenario, hadrons 
(mostly pions), that would have been produced from 
string fragmentation, are converted instead to valence 
quarks and/or antiquarks with current quark masses. In- 
teractions among these partons are again described by 
the ZPC parton cascade model. Since inelastic scatter- 
ings are not included, only quarks and antiquarks from 
melted strings are present in the partonic matter. The 
transition from the partonic matter to the hadronic mat- 
ter is achieved using a simple coalescence model, which 
combines two nearest quark and antiquark into mesons 
and three nearest quarks or antiquarks into baryons or 
anti-baryons that are close to the invariant mass of these 
partons. The present coalescence model is thus some- 
what different from the ones recently used extensively 
[33L I34L I35L l3t| for studying hadron production at in- 
termediate transverse momenta. Using parton scattering 
cross sections of 6-10 mb, the AMPT model with string 
melting was able to reproduce both the centrality and 
transverse momentum (below 2 GeV/c) dependence of 
the elliptic flow ^} and pion interferometry [22j mea- 
sured in Au+Au collisions at y/s =130 AGeV at RHIC 
[37l l38| . It has also been used for studying the kaon in- 
terferometry j3^ and charm flow |2(j in these collisions. 
We note that above parton cross sections are significantly 
smaller than that needed to reproduce the parton ellip- 
tic flow from the hydrodynamic model |4l|. The result- 
ing hadron elliptic flows in the AMPT model with string 
melting are, however, amplified by modelling hadroniza- 
tion via quark coalescence [3(j, leading to a satisfactory 
reproduction of experimental data. As shown earlier 
in Refs. [42L |43|. the hadronization of the quark-gluon 
plasma leads in general to such acceleration at freeze- 
out. 



angle with respect to the reaction plane |4J], i.e. 



III. PSEUDORAPIDITY DEPENDENCE OF 
ANISOTROPIC FLOWS 



The anisotropic flows v n of particles in heavy ion col- 
lisions are the Fourier coefficients in the decomposition 
of their transverse momentum spectra in the azimuthal 
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Because of the symmetry 



(2) 

in the collision geom- 
etry, sine terms do not appear in above expansion. Since 
the projectile and target sides of the reaction plane can in 
principle be identified experimentally, all isotropic flows 
VnS are generally finite. Although in collisions with equal 
mass nuclei the odd-order anisotropic flows of particles 
vanish as a result of the symmetry (y, <p) — ► (— y, (f> + tt). 
at midrapidity and midrapidity only, in case of collisions 
of asymmetric nuclei the odd flow components are finite 
at midrapidity also. 

Anisotropic flows depend on both particle transverse 
momentum and rapidity, and for a given rapidity the 
anisotropic flows at transverse momentum pt can be 
evaluated according to 



Vnipr) = (cos(n4>)} 



(3) 




where (• • •) denotes average over the azimuthal distribu- 
tion of particles with transverse momentum px- 

The anisotropic flows v n can further be expressed in 
terms of single-particle averages: 



(4) 
(5) 

(6) 

(7) 



where p x and p v are, respectively, projections of particle 
momentum in and perpendicular to the reaction plane. 
Since the AMPT model also provides information on the 
spatial anisotropy of colliding nuclear matter, which is re- 
sponsible for generating the momentum anisotropic flows, 
it is of interest to introduce the spatial anisotropic coeffi- 
cient s n by expressions similar to those for the anisotropic 
flows v n but in terms of the spatial distributions of par- 
ticles in the transverse plane. 

In Figs, n (a) and Cfb), we show the pseudorapid- 
ity dependence of v\ and U3, respectively, for charged 
hadrons from minimum bias events of Cu+Au collisions 
at y/s — 200 AGeV by using the string melting sce- 
nario with parton scattering cross sections a p = 3 (open 
squares) and 10 mb (solid squares) and also the scenario 
without string melting (default AMPT model, solid tri- 
angles) . Compared with the AMPT results shown in Fig. 
2 of Ref. |24j for Au+Au collisions at same cm. energy 
per nucleon, where the string melting scenario with a 
parton scattering cross section of 10 mb describes very 



3 




all three scenarios considered here. 
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FIG. 1: (Color online) Pseudorapidity dependence of vi (a) 
and V3 (b) for charged hadrons from minimum bias events 
of Cu+Au collisions at *Js = 200 AGeV by using the string 
melting scenario with parton scattering cross sections a p = 
3 (open squares) and 10 mb (solid squares) and the default 
AMPT model without string melting (solid triangles). 



well the experimental data around midrapidity from the 
STAR collaboration ^(|, asymmetric Cu+Au collisions 
display clearly a stronger v\ around mid-77 in the AMPT 
model with string melting. Also, the pseudorapidity de- 
pendence of v\ has a negative slope around mid-77 and is 
sensitive to the parton scattering cross section. It is fur- 
ther seen from Fig. Q]( a ) that the pT-hxtegrated V\ is very 
small around mid-77 m the default AMPT model without 
string melting. In addition, v\ displays an asymmetry 
in forward and backward rapidities around mid-77 with 
the magnitude of v\ at forward rapidity larger than that 
at backward rapidity. As shown later, this asymmetry 
is related to the different spatial deformation at forward 
and backward rapidities around mid- 77. For the scenario 
without string melting, an even stronger asymmetry of 
Vi in forward and backward rapidities is seen at larger 
pseudorapidity (I77I > 3), with the magnitude of V\ at 
large backward pseudorapidity (Au-like rapidity) about 
a factor of 2 larger than that at large forward pseudora- 
pidity (Cu-like rapidity). This result can be understood 
from the fact that hadronic degrees of freedom dominate 
at large pseudorapidities as the scenario without string 
melting has been shown in Ref.[24( to give a better de- 
scription of the STAR data , and there are thus more 
hadrons and hadronic rescattering at the Au-like rapid- 
ity. We will return to this point later. On the other 
hand, the magnitude of pr-integrated V3 is very small 
(less than 0.3%) in the whole pseudorapidity region for 
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FIG. 2: (Color online) Same as Fig. Qfor V2 and va. 

The pseudorapidity dependence of tj 2 and W4 for 
charged hadrons from minimum bias events of Cu+Au 
collisions at ^/s = 200 AGeV is shown in Fig. [21 by using 
same three scenarios as in Fig. ^ Similar to v\ shown 
in Fig. ^ the elliptic flow vi displays a clear asymmetry 
in forward and backward rapidities around mid-77 as we H 
as at large pseudorapidity. In the case of string melting 
scenario with a parton scattering cross section of 10 mb, 
the value of tj 2 around mid- 77 is about 3%. Compared 
with the V2 value of about 4.5% around mid-77 m Au+Au 
collisions at same cm. energy per nucleon shown in Fig. 
3 of Ref. [24f], the scaling of elliptic flow according to 
the reaction system size as proposed in Ref. l]| is also 
satisfied in asymmetric collisions. As in the case of V3 
shown in Fig. ^ (b) , the magnitude of higher-order px- 
integrated V4 is also very small (less than 0.3%) in the 
whole pseudorapidity region for all three scenarios. 



IV. pr-DEPENDENCE OF ANISOTROPIC 
FLOWS AROUND MIDRAPIDITY 



More detailed information on anisotropic flows can 
be obtained from the differential anisotropic flows, i.e., 
their pt dependence. For hadrons produced around mid- 
77 in heavy ion collisions at RHIC, the partonic degree 
of freedom is expected to be important and the sce- 
nario with string melting describes very well the observed 
anisotropic flows of charged hadrons [24j]. In Fig. |3J 
we show results obtained with the string melting sce- 
nario on the pr-dependence of V\, 1)2, i>3, and V4 for 
charged hadrons at forward mid-77 (0 < 77 < 1.2, left 
panels) and backward mid-77 (—1-2 < 77 < 0, right pan- 
els) from minimum bias events of Cu+Au collisions at 
y/s = 200 AGeV. Both parton scattering cross sections 
Up = 3 (open squares) and 10 mb (solid squares) are used. 
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FIG. 3: (Color online) pr-dependence of vi ((a) and (b)) V2 
((c) and (d)), V3 ((e) and (f)), and V4 ((g) and (h)) for charged 
hadrons at forward mid-77 (0 < 77 < 1.2, left panels) and 
backward mid- 77 (—1.2 < 77 < 0, right panels) from minimum 
bias events of Cu+Au collisions at -^s = 200 AGeV. 



The sensitivity to the parton cross section is clearly seen 
in tj 2 and V4 at both forward and backward mid- 77. We 
note that the elliptic flow exhibits a stronger sensitivity 
to the parton cross section in asymmetric Cu+Au colli- 
sions than in symmetric Au+Au collisions. For example, 
the ratio of vi with a p = 10 mb to that with a p = 3 
mb is about 1.4 in Cu+Au collisions but is reduced to 
about 1.2 in Au+Au collisions as given in Fig. 2 of Rcf. 
|17| . A similar conclusion is obtained for v<± . Sensitivities 
of the pT-dependence of v\ at forward mid-77 and tj 3 at 
backward mid-77 to the parton cross section are also seen. 

Fig. |31 also shows that the pr-dependence of v\ at for- 
ward and backward mid-77 exhibits very different behav- 
iors. At forward mid-77, ui(pr) is non-zero and changes 
from negative to positive values at a balance transverse 
momentum while i>i(pr) at backward mid- 77 is essentially 
zero at low pt and becomes large and positive above 
about 0.8 GeV/c for both parton cross sections. As a 
result, the pT-integrated v\ of charged hadrons around 
mid-77 displays the asymmetry in forward and backward 
rapidities shown in Figs. H( a )- For v 2, its value at back- 
ward mid-77 is slightly larger than at forward mid-77 while 
the value of V4 exhibits an opposite behavior. The dif- 
ference between the p-r-dependence of V3 at forward and 
backward mid-77 is, on the other hand, not so clear due 
to their small values. 

The asymmetry in forward and backward rapidities 
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FIG. 4: (Color online) Time evolutions of si (a), S2 (b), vi 
(c) and V2 (d) of partons at forward mid-77 (0 < 77 < 1.2, solid 
squares) and backward mid-77 (—1-2 < 77 < 0, open squares) 
from Cu+Au collisions at *fs = 200 AGeV and 6 = 6.74 fm 
in the scenario of string melting with a v — 10 mb. 



around mid-77 observed in Fig. [3] can be better under- 
stood from the time evolutions of parton spatial defor- 
mation and anisotropic flows as the latter are transferred 
to those of hadrons when they are formed from quark 
and/or antiquark coalescence. Also, scatterings among 
hadrons, which are included in the AMPT model, do not 
affect much hadron anisotropic flows as a result of small 
spatial anisotropy and low pressure in hadronic matter 
[l0| . In Fig. 01 we show the time evolutions of si, s%, 
Vi and V2 of partons at forward mid-7/ (0 < 77 < 1.2, 
solid squares) and backward mid- 77 (—1.2 < 77 < 0, open 
squares) from Cu+Au collisions at ^fs — 200 AGeV and 
b = 6.74 fm in the scenario of string melting with a p = 10 
mb. We note that the impact parameter b — 6.74 fm cor- 
responds to the same 6/6 max as in Au+Au collisions at 
6 = 8 fm, where 6 max is the sum of the radii of colliding 
nuclei. Similar results for S3, S4, V3 and V4 are shown in 
Fig.l 

It is seen from Fig. 0] (a) that the strength of the 
spatial deformation coefficient si is initially large and 
decreases with time. It reaches a saturation value at 
about 10 fm/c for partons at both forward and backward 
mid-77. The strength of s\ at forward mid-77 is, however, 
larger than that at backward mid-77, leading thus to a 
stronger vi at forward mid-77 as seen in Fig. 0] (c) and 
also in Fig. Q](a). On the other hand, Fig. 0] (b) shows 
that the strength of parton spatial elliptic deformation 
S2 is slightly larger at backward mid-77 than at forward 
mid-77, resulting in a slightly stronger tj 2 at backward 
mid-77 that is shown in Fig. 0] (d) and also in Figs. |21 
and |3| Therefore, the observed asymmetry of V\ and 
i>2 in forward and backward rapidities around mid-77 is 
due to an asymmetry in initial spatial deformation in 
asymmetric heavy ion collisions. 
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(Color online) Similar as Fig. 2] for S3 (a), S4 (b), H3 
w 4 (d). 



From Fig. 03(a), we observe that the S3 displays differ- 
ent time evolutions for partons at forward and backward 
mid-?7 and that its saturated value at backward mid-77 is 
larger than at forward mid-77. As a result, a stronger V3 
at backward mid-77 as shown in Fig. [S] (c) is seen. For 
S4, it also displays different time evolutions for partons 
at forward and backward mid-77. Its saturated value is, 
on the other hand, larger at forward than backward mid- 
77, leading to a stronger at forward mid-77 as shown 
in Fig. |S] (d) and also in Fig. However, because of 
its small value, the asymmetry in forward and backward 
rapidities is not so clear for the V3 of charged hadrons 
around mid-77 as shown in Figs. [2(b) andJSJ 



Pt-DEPENDENCE OF ANISOTROPIC 
FLOWS AT LARGE RAPIDITY 



For hadrons at large rapidities (\tj\ ^ 3) in heavy ion 
collisions at RHIC, the initial dynamics is expected to be 
dominated by the string degree of freedom as they are 
produced later in time when the volume of the system 
is large and the energy density is low. Indeed, it has 
been shown that the default AMPT model in the sce- 
nario without string melting can describe simultaneously 
data on v% and 1*2 at large pseudorapidity ((77! ^ 3) in 
Au+Au collisions at y/s = 200 AGeV 0. Therefore, it 
is of interest to study anisotropic flows at large rapidi- 
ties in asymmetric Cu+Au collisions at RHIC based on 
the AMPT model without string melting. In Fig. [fJl we 
show the predicted pr-dependence of Vi , v%, V3, and i>4 
of charged hadrons at forward (77 = 4, solid squares) and 
backward (77 = —4, open squares) pseudorapidities from 
minimum bias events of Cu+Au collisions at yfs = 200 
AGeV. It is seen that the directed flow v\(pt) shown in 
Fig. (a) exhibits a strong asymmetry in forward and 
backward rapidities, with charged hadrons at backward 
pseudorapidity having a much stronger v\(pt) than that 



CO 

c 
o 

1_ 

T3 
CO 

sz 

TD 
CD 

O) 
l_ 

CD 
_£Z 

o 

■4— 

o 
> 



(a) v , 1 - 

_ = 1 - n ---n s S _ 


(b) v 

: ^- * : 

Am 

B " 


. ^ X 1 

-AMPT (Default) 

. Mini, bias Cu+Au@200 AGeV \ 

— 7-1=4| ^Tl=r4 pi_ 


( C ) V 3 


-(d) v 4 


1 


1.1.1 



0.0 0.5 1.0 1.5 0.0 0.5 

p T (GeV/c) 



1.0 



1.5 



FIG. 6: (Color online) pr-dependence of t;i (a), vi (b), 7J3 
(c) and Vi (d) for charged hadrons at large forward pseudo- 
rapidity (77 = 4, solid squares) and backward pseudorapidity 
(77 = —4, open squares) from minimum bias events of Cu+Au 
collisions at^ = 200 AGeV. 



at forward pseudorapidity. Furthermore, v\ at forward 
pseudorapidity changes from negative to positive values 
at a balance transverse momentum of about 0.8 GeV/c 
while v\ at backward pseudorapidity changes from posi- 
tive to negative values at a balance transverse momentum 
of about 1.2 GeV/c. A balance transverse momentum for 
v\ has also been observed in symmetric Au+Au collisions 
as shown in Fig. 4 of Ref.[2^|, and its existence has been 
attributed to the presence of transverse radial expansion 
0- 

A clear asymmetry of the elliptic flow v% in forward 
and backward rapidities is also seen in Fig. (b), with 
charged hadrons at backward pseudorapidity showing a 
stronger V2(j>t) than that at forward pseudorapidity. A 
similar asymmetry in forward and backward rapidities 
is observed in the higher-order anisotropic flow v^{pt) 
at higher pt as shown in Fig. (c). The higher-order 
anisotropic flow v^{pt) shown in Fig. El(d) is essentially 
zero and it is not clear if it has an asymmetry in forward 
and backward rapidities. 

The observed strong asymmetry of charged hadron 
directed flow v\ and elliptic flow v% in large forward 
and backward pseudorapidities may be related to the 
charged hadron pseudorapidity density, which is shown 
in Fig. \7\ for minimum bias events of Cu+Au collisions 
at y/s = 200 AGeV from the AMPT model in the sce- 
nario without string melting. It is seen that the charged 
hadron pseudorapidity density is asymmetric in forward 
and backward rapidities. Its value at backward pseu- 
dorapidity (Au-like rapidity) is larger than that at for- 
ward pseudorapidity (Cu-like rapidity), leading then to 
more hadronic rescattering and thus stronger v\ and vi 
at backward pseudorapidity as observed above. 
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FIG. 7: (Color online) Charged hadron pseudorapidity den- 
sity for minimum bias events of Cu+Au collisions at y/s = 200 
AGeV from the AMPT model in the scenario without string 
melting. 



around midrapidity in asymmetric Cu+Au collisions dis- 
play a stronger directed flow v\ and their elliptic flow 
is also more sensitive to the parton cross section used 
in the parton cascade. Furthermore, while higher-order 
t>3 and i;4 are small at all rapidities, both v\ and Vi in 
these collisions are appreciable and show an asymmetry 
in the forward and backward rapidities. This asymmetry 
is present even around midrapidity for v\ . Experimental 
verification of these predictions for asymmetric heavy-ion 
collisions at RHIC will be very useful in understanding 
the dynamics of asymmetric collisions. In addition, ex- 
perimentally identifing the projectile and target sides of 
the reaction plane in symmetric collisions is rather dif- 
ficult. This has led to a considerable delay in identify- 
ing the odd flow components, and the accuracy of the 
measurements of these harmonics is still behind the ac- 
curacy and detail achieved in the studies of the elliptic 
flow. Study of collisions with asymmetry systems can 
contribute to a significant improvement of the directed 
flow studies. 
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